Curvature and isocurvature perturbations in a three-fluid model of curvaton decay 
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We study the evolution of the cosmological perturbations after inflation in curvaton models where 
the non-relativistic curvaton decays into both radiation and a cold dark matter component. We cal- 
culate the primordial curvature and correlated isocurvature perturbations inherited by the radiation 
and cold dark matter after the curvaton has decayed. We give the transfer coefficient in terms of 
the initial curvaton density relative to the curvaton decay rate. 
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I. INTRODUCTION 



Observational results obtained by a variety of recent cosmic microwave background and large scale structure ex- 
periments have greatly improved our picture of the Universe on the largest observable scales. These observations are 
consistent with a primordial spectrum of adiabatic, Gaussian, scale-independent perturbations 

An attractive explanation is that these primordial perturbations in the early radiation-dominated era originate from 
quantum fluctuations in an earlier inflationary era where small scale vacuum fluctuations can be stretched up to 
arbitrarily large scales jlfj. Although it is often assumed that primordial perturbations originate from adiabatic 
fluctuations of a slowly-rolling scalar field driving the inflationary expansion (the inflaton), any light scalar field (with 
effective mass less than the Hubble scale) will acquire an almost scale-invariant spectrum during an almost exponen- 
tial expansion. Whether or not these isocurvature fluctuations in non-inflaton fields produce a primordial curvature 
perturbation depends upon whether or not the local value of that scalar field subsequently affects the local radiation 
density [H|l3. 

In the curvaton scenario |lall4Lll5| a light, weakly coupled scalar field remains decoupled from the inflaton during 
inflation and remains decoupled after the end of inflation. Once the Hubble scale drops below the effective mass 
of the curvaton, the field oscillates and eventually decays. But if the decay is late enough then the density of the 
non-relativistic curvaton may have grown large enough to contribute significantly (more than about 1%) to the total 
energy density of the universe. Then the local radiation density after curvaton decay depends upon the local value of 
the curvaton at the end of inflation and it is the initial isocurvature perturbation of the curvaton field that determines 
the primordial density perturbation. This scenario leads to new possibilities for building particle physics models of 
inflation [f6| and qualitatively different constraints upon the inflationary parameters |17| . 

In the conventional slow-roll inflaton scenario there is a remarkable simplification when calculating observational 
constraints: one can calculate the primordial perturbations directly in terms of the inflaton field fluctuations during in- 
flation, independently of the detailed physics of the end of inflation, reheating and the subsequent radiation-dominated 
era. This is because the inflaton fluctuations describe adiabatic perturbations about the background homogeneous 
cosmology on large scales and one can define a conserved curvature perturbation £ directly from energy conservation 

Em 

By contrast, fluctuations in the curvaton field describe non-adiabatic perturbations about the background cosmology 
and the total curvature perturbation £ is no longer conserved on large scales. The non-adiabatic nature of the initial 
perturbations leaves o pen the possibility of residual isocurvature perturbations and detectable non-Gaussianity in the 
curvaton scenario [13.1151 120L l2lj . One can define individual curvature perturbations for the radiation and curvaton 
fluids that remain approximately conserved so long as the fluids can be assumed to be non-interacting 00. This 
has been used to derive an approximate calculation of the resulting primordial perturbations in the curvaton scenario 
in terms of the curvaton density at the decay time in what has been called the sudden-decay approximation 

However an accurate calculation of the resulting primordial perturbations requires numerical solution of the pertur- 
bation equations. This was performed in j^] for a two-fluid model where non-relativistic curvaton particles decay into 
radiation. In this paper we extend this analysis to a three-fluid system to include the possibility that the curvaton also 
decays into non-relativistic cold dark matter (CDM). This allows us to calculate the residual isocurvature perturbation 
that may be inherited by the matter in this scenario. See Ref. |lq for a numerical solution of the perturbed curvaton 
equations in a field theory description. Analytic estimates of the amplitude of the residual isocurvature perturbation 
in other models for the CDM were given in Ref. [2l| . 

In section [H] we set up the coupled evolution equations for a system of dimensionless variables describing the 
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evolution of the background homogeneous cosmology. We then give the coupled evolution equations for linear density 
perturbations in each fluid on uniform-curvature hypersurfaces in section UTTl In an appendix we present equations 
describing the evolution of curvature and isocurvature perturbations in manifestly gauge-invariant form and discuss 
the appearance of singularities for some gauge-invariant quantities in a system such as this that includes energy 
transfer between fluids. This is seen to be due to a singular choice of hypersurface rather than a breakdown of 
perturbation theory. 

The primordial perturbations that result from an initial curvaton perturbation can be characterised by a transfer 
coefficient that is a unique function of a single dimensionless parameter that specifies the initial density of the curvaton 
relative to the curvaton decay rate. We give an accurate analytic approximation (to better than 1%) for the curvature 
transfer coefficient. The same coefficient also determines non-Gaussianity of the primordial perturbations |l4j and, in 
our three-fluid model, the amplitude of the residual isocurvature perturbation. We are able to derive, as an exact result, 
a simple consistency relation between the ratio of isocurvature to curvature perturbations and the non-Gaussianity 
of the primordial perturbations, which was previously derived only using the sudden-decay approximation. 

We present our conclusions in section Hvl 



II. BACKGROUND EQUATIONS AND DECAY 



In this section we will first review the governing equations of the background quantities in a Friedman-Robertson- 
Walker (FKW) universe and then find the evolution equations for the perturbations for three interacting fluids, 
following Ref. |22j. We begin the analysis after the end of inflation when the vacuum energy driving inflation has 
decayed into radiation. We assume that at that time the curvaton density is much smaller than the total density, 
consistent with the assumption that curvaton fluctuations are initially isocurvature perturbations. 

The radiation (7) is a relativistic fluid with pressure P 1 — p-y/3. We model the other two components, the curvaton 
(a) and the cold dark matter (in), as non-relativistic fluids, P„ = P m = 0. This is a good description of the classical 
field dynamics once the Hubble rate has dropped below the effective mass of the curvaton so that the curvaton field 
begins coherent oscillations. 

The evolution of a flat background FRW universe is governed by the Friedmann constraint and the continuity 
equation, respectively, 

H 2 = —P, (1) 
p = -3H(p + P) , (2) 

where the dot denotes a derivative with respect to coordinate time t, H is the Hubble parameter, and p and P are 
the total energy density and the total pressure, related to the density and pressure of the component fluids by 

^p a =p, ^P« = P. (3) 

a a 

The continuity equation for each individual fluid in the background can be written as p3| 

p a = -3H{ Pa + P a ) + Q a , (4) 

where Q a describes the energy transfer per unit time to the a-fluid. 

We assume the curvaton can decay into both radiation and cold dark matter, described by the two decay rates, Fx 
and r 2 respectively. For simplicity we assume the decay rates are fixed constants. This gives the background energy 
transfer as 

Q* = -(Ti+T 2 ) P<7 , (5) 
Q 1 = F \p & , (6) 

Qm = ^2Pa ■ (J) 

The cold dark matter in this model is assumed to be a non-relativistic product of curvaton decay, which does not 
interact with radiation even at these early times. r2/Ti determines the energy density of cold dark matter particles 
relative to the radiation produced by curvaton decay. For the CDM to remain subdominant until matter-radiation 
equality at T eq ~ leV we require 
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decay ^decay 



Toq <icr 6 . (8) 
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The final inequality follows from requiring that the curvaton decays before primordial nucleosynthesis (Td 0C ay > 
lMeV), after which time any significant entropy production would spoil theoretical predictions of the primordial 
abundance of light elements. 

The background energy transfer must always obey the constraint 

^Q Q = 0. (9) 

a 

It is useful to rewrite the background equations (0J), in terms of dimensionless density parameters 



n„ = ^, ^ = ^, n m = ^, (io) 



which gives the evolution equations 



H 

r 2 



n 7 = n^ + n^n^ - 1) , (12) 



f4 - + (13) 

where a dash denotes the differentiation with respect to the number of e-foldings, N = J Hdt. 
The Friedmann constraint Q then yields 

o CT + n 7 + n m = 1. (14) 

The evolution of the inverse Hubble parameter is governed by 



= l + *«7 UiF • ( 15 ) 



Hence we have an autonomous dynamical system consisting of the four Eqs. I|11H13I) . and 115|) . but the constraint 
(1141) makes one of them redundant. 



A. Numerical solutions 



In Fig^we show a phase-plane which plots the evolution of the curvaton energy density against the dimensionless 
decay rate, < ri/(ri + H) < 1. Note that for T2 <C T_\ the evolution of fl a is effectively independent of T2 and the 
phase-plane shown is identical to that studied in R.ef . [22| . 

All trajectories in the phase-plane in FigO] start close to the origin with SX^^C 1 and H ^> T\. From the limiting 



form of Eqs. i|ll|) and (|14fl for fi 7 ~ 1 we can see that fl a oc H 1 / 2 oc e 
identified by different values of the dimensionless parameter 



Thus different trajectories can be 



- initial 



If we consider the initial conditions for the fluid model to be set when the curvaton begins to oscillate (H u 
then we can write 



3M|j 



T7 



1/2 



(16) 



m„ 



(17) 



The curvaton vacuum expectation value is bounded, (c 2 )i n < Mpj, as we work in the regime where there is no phase of 
curvaton-driven inflation [l7|. Nonetheless p can be large if the curvaton decay rate is slow, e.g., m a /T\ ~ (ilfpi/m CT ) 2 
for gravitational strength decay. 

We show in Figures 2 to 4 examples of the time-evolution of the three fluid energy densities along different trajec- 
tories. Figure Ufa) corresponds to trajectory A of the phase-plane Figure ^ The system begins radiation dominated 
but as time passes the curvaton density, with a dust equation of state, grows relative to the radiation density. The 
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FIG. 1: The phase plot shows the evolution of the dimensionless energy density, Q a , for different values of the dimensionless 
parameter, p, defined in Eq. 1)16^ . varying from 3.15 x 10 3 (top line) to 3.15 x 10 -3 (bottom line). For all of the lines T2 < Ti. 
The evolution of the densities of all three fluids are plotted below in Figs. [5] to 2] for the trajectories, marked A, B and C. 



energy density of the curvaton then rapidly decays when H ^ T\. In this example there is a period of curvaton 
domination ~ 1) before the curvaton decays. 

Figure Efa), corresponding to trajectory B of the phase-plane figure QJ, has a lower value the parameter p, so 
the curvaton density does not climb as high before it decays away. Figure Ufa), corresponding to trajectory C of the 
phase-plane figure 0J, has the lowest value of p of the three plots and the ratio of initial curvaton density to the 
decay rate is low enough that there is no period of curvaton domination in this example. 



III. PERTURBED EQUATIONS 



The evolution equation for density perturbations in the a fluid is given on large scales and on uniform-curvature 
hypersurfaces by | 22| 

S Pa +3H(S Pa + SP a ) + Q a ^-SQ a = 0, (18) 

2p 

where 5p a , 5P a and 8Q a are the density perturbations, the pressure perturbation, and the energy transfer perturbation 
of the a-fluid. Sp is the total density perturbation on the uniform-curvature hypersurfaces and appears in the equation 
for the perturbed energy transfer due to the gravitational time dilation . 

The perturbed energy transfer (assuming the physical decay rates are unperturbed) is given by, 

5Qa = -(ri + r 2 )<fAx, (is) 

<5Q 7 = T X 5 P<J , (20) 

SQm = r 2 s Pa . (2i) 



We will evolve the perturbed conservation equations (|18fl . which can be written in terms of the dimensionless density 
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FIG. 2: (a) The evolution of f2 CT (solid), fi 7 (dashed) and fi m (dotted) against e-foldings N of line A of Fig0 The initial 
value of f7 CT is 10" 2 ' 5 , and IT = 10" 10 and T 2 = 1CT 12 . (b) The evolution of Cr/C^in (dotted), in (dot-dashed), Cm/Ciu 

(dashed) and the total curvature perturbation, C/Cr,in (solid), against e-foldings N of line A of FigtTI 



parameters as 



^ + (3 + ^)^-^O CT |=0, (22) 
6^ + 4S Pl - ^ (s Pa - nj-fj = , (23) 
S/4 + 3fy m - ^ fjyo CT - = . (24) 



H \ 2 , 

where a prime denotes derivatives with respect to iV. 

Written in this form it is evident that the linear evolution of the density perturbations depends only upon the 
dimensionlcss background variables Q a and decay rates T a /H. Hence the final amplitude of the density perturbations, 
5 Pa, final, relative to the initial curvaton perturbation <Sp ff ,mitial will depend only upon which trajectory is followed in 
the background phase-space (Q a ,T a / H). They will not depend on dimensional quantities such as -ffmitiai- In fact 
we will show that the large-scale perturbations after the curvaton has decayed can be written simply in terms of the 
initial curvaton perturbation and a function of the dimensionless parameter p defined in Eq. (|16|) . 
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FIG. 3: (a) The evolution of Q a (solid), 57 7 (dashed) and Q m (dotted) against e- foldings N of line B of Fig0 The initial value 
of fi CT is lO -4 ' 6 . Ti and T 2 are the same as in Fig|5] (b) The evolution of CoVCo-in (dotted), C?/ Co .in (dot— dashed), Cm/Co-, in 
(dashed) and the total curvature perturbation, C/Cr,in (solid), against e-foldings N of line B of Fig0 



A. Curvature and isocurvature perturbations 



The gauge- invariant variable, £, introduced by Bardeen |24l l25j . is commonly used to describe the large-scale 
curvature perturbation. It is equivalent to the curvature perturbation of uniform-density hypersurfaces and remains 
constant on super-horizon scales for adiabatic perturbations. This quantity is simply related to the total density 
perturbation on uniform-curvature hypersurfaces flil l25| 
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(25) 



We can also define the curvature perturbation on uniform a-fluid density hypersurfaces, £ Q , which is related to the 
a-fluid density on uniform-curvature hypersurfaces [Isf 



Pa 



(26) 



The total density perturbation is just the sum of the density perturbations of the individual fluids, 5p = ^ Q Sp a , and 
hence the total curvature perturbation £ is the weighted sum of the individual perturbations [18( , 
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FIG. 4: (a) The evolution of Q a (solid), 57 7 (dashed) and Q m (dotted) against e- foldings N of line C of Fig0 The initial value 
of n a is 10" 5 ' 5 . Ti and T 2 are the same as in Fig|5] (b) The evolution of Co-/Cr,m (dotted), Ct/Cct-.u! (dot— dashed), Cm/Co-, in 
(dashed) and the total curvature perturbation, C/Cr,in (solid), against e-foldings N of line C of FisTTl 



The difference between any two individual curvature perturbations, £ a and (^3, describes an isocurvature or entropy 
perturbation OHl 

S a p = 3 (C« - C/s) ■ (28) 

Note that the definition of Q a becomes singular whenever p a — 0. For instance the physical density of cold dark 
matter is initially vanishing, but growing due to the decay of the curvaton (p m > 0). However after the curvaton has 
completely decayed, the matter density is diluted by the expansion (p m < 0). At one instant we have p m = and 
C m is ill-defined. This has been previously noted in the case of an oscillating scalar field at stationary values of the 
field p^ |. It is a consequence of the uniform a-density hypersurfaces becoming ill-defined, and does not signal the 
breakdown of perturbation theory, as long as we can find some gauge where the hypersurfaces are well-defined and 
perturbations remain small (see Appendix). 

In this case we will numerically solve for the fluid density perturbations on uniform-curvature hypersurfaces. Fig. 
shows an example of the evolution of the density perturbations in equations lf2*2^) to ij2*4*|) above, showing that they 
do indeed remain well-defined throughout the evolution. 

From the density perturbations we then construct the curvature perturbations, Eci. (|25[) and (|26|) . which gives 
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FIG. 5: The evolution of the density perturbations on uniform-curvature hypersurfaces, b~p a (solid), 5p 7 (dashed) and 8p n 
(dotted), against e- foldings N along the background trajectory B in Fig0 
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At late times, after the curvaton has completely decayed the fluids are non-interacting and have monotonically 
decreasing energy densities (diluted by the expansion) so the curvature perturbations, £ a , are well-defined, and 
constant on large scales. Observations of, for example, the CMB directly constrain the amplitude of the total 
curvature and isocurvature perturbations. On large an gula r scales the temperature perturbations on the surface of 
last-scattering are given by the Sachs- Wolfe effect 0, 0, |2Z| 
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(32) 
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B. Numerical results 



We have numerically evolved the coupled system of background equations and 115fl and first-order per- 

turbations H22fl -124 [) . In each case we begin with a smooth (C 7 ,in = 0), dominant radiation fluid (CI, 



7i ln 



1) and 

a subdominant curvaton fluid (f2 CTj i n <C 1) with an initial perturbation Co-, in ^ 0. The cold dark matter density is 
initially zero (fl m ^ n = 0, Cm, in = 0), and is produced solely from curvaton decay. The initial density perturbations on 
uniform-curvature hypersurfaces are given in terms of Ca.in from Eqs. I|29|) - (|31|l . 

In all cases we start the evolution at some high energy scale, with a large value of the Hubble parameter, characterised 
by a small initial value of the Hubble time relative to the decay time, Ti/Hi n <C 1. The subsequent evolution is not 
directly affected by the precise value of ri/iJ; n , but rather by the parameter p defined in Eq. l|lt)fl which characterises 
the phase-space trajectory. In particular the final values of the curvature and isocurvature perturbations are a function 
solely of the parameter p. 

Figures Hfb), E^b) and^Jb) show the evolution of the large-scale curvature perturbations for the three fluids along 
three different background trajectories. 



Radiation 



The radiation curvature perturbation C 7 grows as the curvaton decays, and reaches an equilibrium point in the 
phase-space where £ 7 = Co- 22]. The final value for £ 7 relative to the initial value of Co-, in depends upon the relative 
density of the curvaton at the decay epoch, f^dccay This determines the proportion of the final radiation that comes 
from the decay of the perturbed curvaton. 
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FIG. 6: Curvature transfer coefficient r = C7, out /Co-. in, plotted against background parameter p = jfi CT (H/Ti) 1 ^ 2 



For example, in Fig.^a), corresponding to trajectory A of Fig. ^ the curvaton comes to dominate (£l a ~ 1) before 
the decay epoch. Thus almost all the radiation at late times comes from the decay of the curvaton and hence, in 
Fig. |2Ib), the radiation fluid at late times inherits the initial curvaton perturbation, £ 7 ~ Cr.m- Compare this with 
Fig. El corresponding to trajectory C of the Fig.^ In this case the curvaton decays while ft a <C 1 and the proportion 
of the final radiation energy density which results from the decay of the perturbed curvaton component is much 
smaller in this picture. 

We define r to be the ratio between the curvature perturbation of the radiation component at the end of the 
calculation (when all of the curvaton component has decayed away) to the initial curvaton perturbation, 



c 



7, out 



We find that r is solely a function of the choice of trajectory in the background phase-space, Fig. QJ. 
is uniquely identified in terms of the initial data by the value of the parameter p defined in Eq. (|16|l 
We plot r against p in Fig. ©. At small values of p this approaches a linear function 

r ~ 0.924p , 



(33) 

Each trajectory 
(34) 



see Fig. {7J). This coincides with the sudden-decay approximation calculated in Ref. [22] for fl a _i n <C 1. The linear 
approximation is good to better than 1% up to p = 0.011 (r = 0.0108) and better than 10% up to p = 0.11 (r = 0.10). 
At larger values of p the transfer coefficient r{p) turns over and r — > 1 for p> 1. The full function is well fit by 



ffit(p) = 1 



0.924 > 9 
1 H p 



(35) 



where, by fitting over the range 10 2 < p < 10 2 , we find q 
to better than 0.2% for any value of p. 



1.24. This gives a fit to the numerical calculation of r(p) 



2. Cold Dark Matter 



In our model the universe starts with an initial cold dark matter density, p m , equal to zero which then grows 
due to curvaton decay, but eventually starts decreasing due to cosmic expansion. As remarked earlier the curvature 
perturbation associated with the cold dark matter, £ m , becomes ill-defined when p m = 0. This occurs when Q rn = 
T2Pa — 3i?p m , in Eq. 10} for a — m. The divergence of £ m is clearly seen at early times (before the curvaton decays) 
in Figures |1>) , [|b) and H» . 

The curvature perturbation £ m is initially time-dependent due to the energy transfer from the perturbed curvaton, 
but settles down to a constant value at late times, once the curvaton has completely decayed ri + T2 ^> H . We find 
that if the cold dark matter is produced solely from curvaton decay, it inherits a curvature perturbation Cm, out at late 
times equal to the initial curvaton perturbation, Cer.in- 




FIG. 7: Curvature transfer coefficient r, defined in Eq. 13311 . plotted against parameter p, defined in Eq. If 6H . with a linear fit 
r = 0.924p for p <C f . This linear fit is accurate to within f % for p < O.Of f . 




FIG. 8: Isocurvature transfer coefficient, 1 — r, defined in Eq. (I33H . plotted against parameter p, defined in Eq. If 6H . Also shown 
(dot-dashed line) is the fitting function f — rat, defined in Eq. I|35|l . This is accurate to within 1% for p < 8.14 (corresponding 
to I - r > 0.086). 



This simple result can be derived by considering a composite energy density constructed from the CDM density 
and a fixed fraction of the curvaton density: 

r 2 

Pcomp = Pm + z~ —fTPa ■ (36) 

f 1 + 1 2 

From Eqs. © and JJJ) we see that the corresponding energy transfer is zero, i.e., 

r 2 

Qcomp = Qm "i" ~ ' ~^~Qa — • (37) 
1 1 + 1 2 

Thus we obtain from Eqs. I0J and l|18l) the standard evolution equations for the composite energy density and its 
perturbation: 

pcomp — -SHp comp ) (38) 

^Pcomp = -3H6p comp . (39) 

Note that as both the curvaton and CDM are pressureless then so is the composite energy density, Pcomp = 0. Since 
we have a conserved energy density with a unique equation of state the corresponding curvature perturbation, defined 
by Eq. ®, 



C - ^ com P ( |.m 

*3Pcomp 
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is conserved on large scales . 

Initially p m = and Sp m = and hence Ccomp = ( aM . At late times p m > (r 2 /(ri + T 2 )) p a , and Ccomp = Cm ,out- 
Therefore 

Cm ,out = Co",m J (41) 

as we see in Figs. |2|3]and0] 



3. Curvature and isocurvature perturbations 

After the curvaton has decayed, (fi CT — ► 0), we are left with one overall curvature perturbation, Eq. ((27(1 . 

C 4 ^ C I 3 ° m C (42) 

^ ~ 4fi 7 + 3fi m ^ 7 + 4ft 7 + 3fi m U ' 1 j 

and one relative isocurvature (entropy) perturbation between the radiation and cold dark matter, Eq. 1)28(1 . 

5 m7 = 3 (C m - Cr) • (43) 

The primordial curvature perturbation is usually taken to be the curvature perturbation at, for instance, the epoch 
of primordial nucleosynthesis when the Universe is known to be radiation dominated (f2 7 3> f2 m , tt a — 0), and hence 

C^C 7 - 

We have seen that the primordial curvature perturbation relative to the initial curvaton perturbation is a function 
of the trajectory, r(j>) in Eq. 1(33(1 . 

Clnuc = Cy,out — r Cff,in , (44) 

whereas the final cold dark matter curvature perturbation is identical to the initial curvaton perturbation, Eq. 141|) . 
Hence 

S \nuc = Sm 7 ,out = 3 (l - r) C<r,in ■ (45) 

If r = 1 then we obtain an adiabatic curvature perturbation after the curvaton has decayed and there is no late-time 
isocurvature perturbation. 

Note that the fitting function ratip) given in Eq. 1)35(1 gives a good approximation, for 1 — r (better than 1% for 
r < 0.914 (p < 8.14) or better than 10% for r < 0.993 (p < 66)). 

Therefore the primordial isocurvature and curvature perturbations are 100% correlated with their ratio fixed by 
r{p): 



(46) 



This expression was derived in Ref. 20] using the sudden-decay approximation. We have shown that it is in fact an 
exact expression when the CDM is produced directly from the decay of the curvaton, and r is defined by Eq. 1(33(1 . 



IV. CONCLUSIONS 



In this paper we have studied the coupled evolution of background densities and linear perturbations of three 
interacting fluid components: radiation, a curvaton field and cold dark matter. We solve numerically the evolu- 
tion equations for density perturbations on spatially flat hypersurfaces to calculate the spectrum of curvature and 
isocurvature (or relative entropy) perturbations that result from initial perturbations in the curvaton field. 

Whereas the resulting radiation curvature perturbation, Cr.out, is dependent on the initial values of background 
parameters (such as curvaton density), we find the simple result that the final CDM curvature perturbation, Cm. out, 
is always exactly equal to the initial curvaton perturbation Cer.in, in a model where the CDM is directly produced by 
curvaton decay. We are able to demonstrate this analytically by constructing a conserved quantity, Ccomp, that at 
early times is equal to the curvaton perturbation and at late times is equal to the CDM curvature perturbation. In 
effect this solves the evolution of one degree of freedom and allows us to solve numerically only for C 7 as a function 
of input parameters. 
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The background evolution of the fluids is described by a trajectories in phase-space distinguished by a dimensionless 
parameter p = ^(H/Ti) 1 / 2 }^. The ratio r between ( 7lOU t and Ccr,in is a function solely of this trajectory parameter p 
and is well-fit by a simple numerical function given in Eq. H35f) . This coincides with the "sudden-decay approximation" 
00 forp< 1. 

The resulting radiation and CDM perturbations after curvaton decay can then be used to predict the spectrum 
of CMB anisotropics in a given model where the CDM is produced directly as a curvaton decay product. CMB 
constraints are usually quoted in terms of the primordial curvature £ nuc and an isocurvature perturbation 5 nuc around 
the epoch of primordial nucleosynthesis, leading to 



(47) 



In the curvaton scenario (in = C7,in = and S- m — 3(( CT) i n — (7, in) = 3( CT .j n . Although the model we have studied 
describes the production of CDM directly from the curvaton decay, one would expect a similar result to hold for 
models where the net baryon or lepton number is directly produced from curvaton decay. 

The matrix equation (|47|) is a particular case of the general form of the transfer matrix for curvature and isocurvature 
perturbations presented in Ref. |(J |28| where the contribution of the initial isocurvature perturbation to the resulting 
primordial curvature and isocurvature perturbation is model dependent. For instance in a curvaton scenario where 
the CDM is a thermal relic which decouples from radiation sometime after the curvaton has decayed we have instead 



s) -at)( ( s 



(48) 



The matrix equations H47[) and (|48|1 using the numerically determined function r(p) [or its approximate form r$x(p)] 
enable one to relate CMB observations to curvaton model parameters such as the curvaton decay rate, T, and the 
curvaton density, at some early time, characterised by the value of the Hubble rate H, and the initial spectrum 
of curvaton perturbations ( CT = (l/3)Sp (T /p cr . 

For example, in the model where the CDM is directly produced by curvaton decay we confirm [T^j that the relative 
amplitude of primordial isocurvature to curvature perturbations is a fixed ratio 3(1 — r)/r which is a function solely 
of the background parameter p defined in Eq. i|16f) . This holds independently of the "sudden-decay approximation" 
previously used. 

In any curvaton model the transfer coefficient r is independent of scale in the large-scale limit (by definition) so 
that the scale-dependence of the primordial perturbation spectra are inherited directly from the scale dependence of 
the curvaton perturbation at the end of inflation 0, |2S| 

An CT = - ° a = -2e + 2r h . 49 
a In k 

In this expression e = —H/H 2 is the usual dimensionless slow-roll parameter during inflation and r\ a = m 2 /3H 2 is 
the effective mass of the curvaton relative to the Hubble scale during inflation. Both the curvature and isocurvature 
perturbation (if any) originate from initial curvaton perturbation and are thus completely correlated, independent of 
the physics of the curvaton decay. In the absence of any additional source for the isocurvature perturbations, the 
curvature and any isocurvature perturbations must share the same spectral tilt, given by Eq. (|49|l . 

Note that the curvature transfer coefficient r can also be related to the non-Gaussianity of the primordial perturba- 
tions. If r is small then the initial curvaton perturbation Ca,in must be correspondingly large in order to produce a given 
amplitude of primordial curvature perturbation. This leads to a larger contribution from second-order perturbations 
in the curvaton density, which shows up in the non-linearity parameter 0, |29| 

/nl^, (50) 

for r < 1 0| . (See |3(j for corrections when r ~ 1.) There is thus a possible consistency test of this curvaton 
model using two observables. Existing constraints on the relative amplitude of primordial isocurvature to curvature 
perturbations @ constrain r > 0.9 in models where the CDM is directly produced by curvaton decay, and hence limit 
the allowed non-Gaussianity. However in curvaton models which produce purely adiabatic perturbations (<S nuc = 0), 
determinations of the non-linear parameter /nl may be the only way to determine independently both the transfer 
parameter r (and hence the model parameter p) and the initial amplitude of the curvaton perturbation Q a . 
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APPENDIX A: GAUGE-INVARIANT ADIABATIC AND ENTROPY PERTURBATIONS 

In a multi-component system it is often useful to re- write the coupled perturbation equations in terms of the instan- 
taneous overall curvature perturbation (often called the adiabatic perturbation) and the relative entropy perturbations 
between fluids. In this appendix we follow the treatment of Ref. [221 ] and the reader is referred to this source for more 
details. Throughout this paper we work in the large scale limit, i.e., neglecting gradient terms. 
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1. Evolution of curvature and entropy perturbations 

The gauge- invariant definition of the total curvature perturbation on uniform-density hypersurfaces is 0, l24l l3l| 

C = -^-H S 4, (Al) 
P 

where ip and Sp are the gauge-dependent curvature and density perturbations. The evolution of the curvature 
perturbation is given on large scales by 0, I22I l32| 

C = ^T^nad , (A2) 

p + P 

where the non-adiabatic pressure perturbation is <5P na d = SP — cj;Sp and the adiabatic sound speed is Cg = P / p. Thus 
the total curvature perturbation is constant on large scales for purely adiabatic perturbations. 

In the presence of more than one fluid, the total non-adiabatic pressure perturbation, <5P na d , may be split into two 
parts, 

5P aad = 5P int r + 5P rcl . (A3) 
The first part is due to the intrinsic entropy perturbation of each fluid 

<5Pi„tr = ^<5Pintr,a, (A4) 

OL 

where the intrinsic non-adiabatic pressure perturbation of the a-fluid is given by 

<5Pntr,a = 5P a - (? a &p a , (A5) 

where = P a /p a is the adiabatic sound speed of that fluid. The second part of the total non-adiabatic pressure 
perturbation i|A3|l is due to the relative entropy perturbation between different fluids, 

where S a is the relative entropy (or isocurvature) perturbation, Eq. Q28JI . 

Sap = -3H( S 4Z- 6 -^). (AT) 

V Pol P(3 J 

Note, that Eq. l|A6f) corrects a sign error in Eq. (2.31) of Ref. pa |. 

The evolution equation for the relative entropy perturbation S a p on large scales is given by p^ 

A „ „ / 3H5Pi n tr t0 , — SQi n t r a — (JQrel.a 3if (JPintr.fl — <%3intr,3 — ^Qrel,0 \ / a o \ 

O a 8 = oil ■ ; ■ ; , (AS) 

V Pa PP J 

where the intrinsic non-adiabatic energy transfer perturbation is defined as p^ 

5Qintr,a = SQ a - ^-5p a , (A9) 
Pa 

and the relative non-adiabatic energy transfer is defined as 

xn QaP ( Sp a Sp\ 0« V^. c , A1 „s 

2p V Pa P J ®Hp ^ 

For fluids, such as radiation or non-relativistic matter in our curvaton model, with definite equation of state, 
Pa(pa), the intrinsic non-adiabatic pressure perturbation in Eq. l|A4(l is zero. Hence the total non-adiabatic pressure 
perturbation (IA6|) for our model can be written as 

£Piad = 5P Icl = —^7 (/9 m <S m7 + p a S&~) . (All) 
\)tl p 
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Hence the evolution equation i|A2(l for £ is 



C — ~3^T (pa^ay + PmSmy) ■ 



The intrinsic energy transfer perturbation, following from Eq. (|2ip. are 





= o, 




r x . 




r 2 . 



while the relative energy transfer perturbations are 

(Ti + r 2 )p ff 



£Qrel,er 



6Hp 



6ii> 



(Pa^ay ~\~ PmSm-f*) ; 
{Pa^am P^/Smj*) ■ 



Hence the evolution of the relative entropy perturbations can be written as 



1 

2p" 



1 \Pa—r~ «->ni7 



P- 



(Ti + r 2 ) / o C r^-5 (Tm 

Pa 



Pn 



c — I r p I ^ c r ( i ''" \ '"" c r ( i '"" 1 '''' c 

"mi — 1 I" 1 2~ — TT^nn + 1 2 1 — — - — £> am — 1 1 I 1 — — —0,77 , 



A) 



At \ P. 



2Y x p- T 2 p a ^ L 

Pl Pa 



Po 



-(r 1 + r 2 )^^ 7 

2P Pa 



To 1 



2/5/ Pn 
Pa \ Pa 



2pJ p 



(r x + r, 



jPa_Pn 
2p Pa 



2pJ p\ 



Pa Pj 



2 „ . ^7 ■ 
2/3 p m 



(A12) 

(A13) 
(A14) 

(A15) 

(A16) 
(A17) 
(A18) 

(A19) 
(A20) 
(A21) 



Note that we can always identify a trivial solution to Eqs. l|A19llA"2TJ> with all the entropy perturbations S a p — 
which corresponds to the adiabatic mode and leaves £ = constant in Eq. l|A12fl . even in a system such as this including 
energy transfer between different components. 



2. Singular and Non— Singular Hypersurfaces 

From Eq. (|A7(1 we see that the relative entropy perturbation S a p becomes singular whenever p a — 0. This can 
occur if we have multiple fluids with energy transfer between them whenever Q a — 3H(p a + P a ). This is because 
the uniform a- fluid hypersurfaces become ill-defined to first-order whenever p a is stationary, i.e., the first-derivative 
vanishes. 

Thus a divergence in a gauge-invariant quantity such as S a p need not signal the breakdown of perturbation theory, 
but instead may indicate that a particular choice of hypersurface was not well-defined [2(| . This is indeed the case in 
the curvaton model studied in this paper: we see that the density perturbations on uniform-curvature hypersurfaces 
remain finite even when p a = 0. 

From the definition of <S Q ^, Eq. (|A7|1 above, we see that we can construct a non-singular quantity 

= PaPpSap = -3H (pp5p a - PaSpp) • (A22) 

Since the definition of S a p is gauge invariant then so is A Q ^. 

The evolution equation for the total curvature perturbation on large scales, Eq. (|A2|I can then be written in terms 
of A a j3 as 



1 



3H(p + P) 



1 Z_^i 



H5P intr , a - c 2 a A «/3 \ ■ (A23) 



The uniform total density hypersurfaces, on which £ is defined, always remain well defined, since p ^ after 
inflation. 



